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Abstract. In this paper, we describe and study a mathematical framework 
for cryo-elecron microscopy. The main result, is a a proof of the admissability 
(correctness) and the numerical stability of the intrinsic reconstitution algo- 
rithm which was introduced by Singer and Shkolnisky in [7 . In addition, we 
explain how the various numerical observations reported in that work, follow 
from basic representation theoretic principles. 



0. Introduction 

"Three dimensional electron microscopy" pQ is the name commonly given to 
methods in which the three dimensional structure of a macromolecular complex 
is obtained from the set of images taken by an electron microscope. The most 
general and widespread of this methods is single-particle reconstruction. In this 
method the three dimensional structure is determined from images of randomly 
oriented and positioned identical macromolecular complexes, also referred to as 
molecular particles. A variant of this method is called cryo-electron microscopy (or 
cryo-EM for short) where multitude of molecular particles are rapidly immobilized 
in a thin layer of ice and maintained at liquid nitrogen temperature throughout 
the imaging process. Single particle reconstruction from cryo-electron microscopy 
images is of particular interest, since it promises to be an entirely general technique 
which does not require crystallization or other special preparation stages and is 
beginning to reach sufficient resolution to allow the trace of polypeptide chains 
and the identification of residues in protein molecules O QJ] . 

Over the years, several methods were proposed for single particle reconstruction 
from cryo-EM images. Present methods are based on the "Angular Reconstitution" 
algorithm of Van Heel [5] , which was also developed independently by Vainshtein 
and Goncharov [2]. However, these methods fail with particles that are too small, 
cryo-EM images that are too noisy or at resolutions where the signal-to-noise ratio 
becomes too small. 

0.1. Main results. In [7], a novel algorithm, referred to in this paper as the intrin- 
sic reconstitution algorithm, for single particle reconstruction from cryo-EM images 
was presented. The appealing property of this new algorithm is that it exhibits 
remarkable numerical stability to noise. The admissibility (correctness) and the 
numerical stability of this algorithm were verified in an overwhelming number of 
numerical simulations, albeit, a formal justification was still missing. In this paper, 
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we prove the admissibility and the numerical stability of the intrinsic reconstitu- 
tion algorithm. The proof relies on the study of a certain operator C, of geometric 
origin, referred to as the common lines operator. Specifically, 

• Admissibility, depends, among other things, on the fact that the maximal 
eigenspace of C is three dimensional. 

• Numerical stability, depends on the existence of a spectral gap which sep- 
arates the maximal eigenvalue A max from the rest of the spectrum. 

In this regard, the main technical result of this paper is a complete description 
of the spectral properties of the common lines operator. In the course, we describe 
a formal mathematical framework for cryo-EM which explains how the various 
numerical observations reported in [7] follow from basic representation theoretic 
principles, thus putting that work on firm mathematical grounds. 

The remainder of the introduction is devoted to a detailed description of the 
intrinsic reconstitution algorithm and to the explanation of the main ideas and 
results of this paper. 

0.2. Mathematical model. When modeling the mathematics of cryo-EM its more 
convenient to think of a fixed macromolccular complex which is observed from 
different directions by the electron microscope. In more details, the macromolecular 
complex is modeled by a real valued function on a three dimensional Euclidian 
vector space V ~ M 3 , which describes the electric potential due to the charge 
density in the complex. A viewing direction of the electron microscope is modeled 
by a point x on the unit sphere X — S(V). The interaction of the beam from 
the electron microscope with the complex is modeled by a real valued function R x 
on the orthogonal plane P x = x , given by the Radon transform of (f> along the 
direction x, that is 

Rx («) = / <t>(v + l)dl, 

for every v £ P x , where L x is the line passing thorough x and dl is the Euclidian 
measure on L x . 

The data collected from the experiment is a collection of Radon images R x : 
P x — > R, x £ Xm\ where Ajy C X consists of N points. The main empirical 
assumption is 

Empirical assumption: The points x € Xm are distributed independently 
and uniformly at random. 

We emphasize that, in practice, the embeddings i x : P x > V, x £ Xn are not 
known. What one is given is the collection of Radon images {R x : P x — > R : x 6 A at}, 
where each plane P^, x € Ajv should be considered as an abstract Euclidian plane. 

The main problem of cryo-EM is to reconstruct the orthogonal embeddings 
{ix ■ Px ^ V : x S Xn} from the Radon images {R x : P x — > M : x G Xjy}. We will 
refer to this problem as the cryo-EM reconstruction problem and note that, grant- 
ing its solution, the function (f> can be computed (approximately) using the inverse 
Radon transform. 

0.3. The Fourier slicing property and the common lines datum. The first 
step of the reconstruction is to extract from the Radon images a certain linear 
algebra datum, referred to as the common lines datum, which captures a basic 
relation in three dimensional Euclidian geometry. 
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First one notes that for every x € Xn, the following relation holds 
(0.1) R x <>■, ■ 

where the operation (— ) on the left hand side denotes the Euclidian Fourier trans- 
form on the plane and the operation (— ) on the right hand side denotes the 
Euclidian Fourier transform on V. This relation follows easily from the standard 
properties of the Fourier transform and is sometimes referred to as the Fourier 
slicing property (see [6]). The key observation, first made by Klug (see [4]), is that 
(jO.ip implies, for every pair of different points x,y £ Xn, that the functions R x 
and R y must agree on the line of intersection (common line), that is 

-Rx|p x nP H = Ry\P x nP v - 
Hence, if the function <fi is generic enough then one can compute from each pair 
of images R x and R y the operator Cn (x,y) : Py — * P x; which identifies the line 
of intersection between the two planes. Formally, this operator is given by the 
composition C XiV o Cy X , where C x , y and C VjX are the tautological embeddings 

C x ,y ■ Fx Py c y f*cC5 
Cy^ x ■ Fx ^ Py e * Py - 

0.4. The intrinsic reconstitution algorithm. The intrinsic reconstitution al- 
gorithm reconstructs the orthogonal embeddings {i x : P x — ► V : x £ Xn} from the 
common lines datum {Cn (x,y) : (x,y) € Xn X Xn}. The crucial step is to con- 
struct an intrinsic model of the three dimensional Euclidian vector space V which 
is expressed solely in terms of the common lines datum. The construction proceeds 
in four steps: 

Ambient vector space: We define the 2N dimensional Euclidian vector space 

n N = p.. 

x£X N 

Common lines operator: We define the symmetric operator Cn ■ T~Cn — » 
H N by 

C N (s)(x) = v ±- £ C N (x,y)( S (y)). 

Intrinsic model: We define the Euclidian subspaceQ 

V w = H N (A) , 

A>l/3 

where Hn (A) denote the eigenspace of Cn associated with the eigenvalue 
A. 

Intrinsic maps: For every x G Xn, we define the map 

¥> x = V / 273-(Pr ie ) t :P !1I -V JV , 
where Pr^ : Yn Px is the orthogonal projection on the component P^. 

The fact that the vector space Vjv is of the right dimension is granted by the 
following theorem: 



^The condition A > 1/3 in the definition of Vjv will be clarified when we will discuss the 
spectral gap property in the next subsection. 
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Theorem 1. For sufficiently large N , we have 

dimVjv = 3. 

The fact that the collection of maps {ip x : P x — * Vn ■ x G Xn} solves the cryo- 
EM reconstruction problem is the content of the following theorem: 

Theorem 2. There exists an (approximated) isometry tn : V « Yn which satisfies 
the following property: 

T N oi x = tp x , 

for every x 6 Xn- 



Remark 1. Theorem^ implies that the vector space V equipped with the tautological 
embeddings {i x : x G Xn} is isomorphic to the intrinsic vector space Yn equipped 
with the mappings {tp x : x G Xn}- Hence, for all practical purposes they are indis- 
tinguishable. This is, to our judgement, an elegant formal example, realizing the 
general phylosophy about the appearence of strucuture form data. 

0.5. Analytic set-up. The proofs of Theorems Q] and O are based on an approxi- 
mation argument of the discrete from the continuos, which we are going to explain 
next. 

Let — > X be the vector bundle on the unit sphere whose fiber at a point x G X 
is the plane = x and let H denote the vector space of smooth global sections 
r (X, f)). The vector space TL is equipped with a inner product, given by 



(si,s 2 )= J B (si (x) , S2 (x)) dx, 

where dx is the Haar measure on the unit sphere. The common lines datum can be 
used to form a kernel of a symmetric integral operator C : H — > H which is given 

by 

C(s)(x)= J C{x,y){s{y))dy, 

The main difference from the discrete scenario is that, here, the space Ti. supports 
a representation of the orthogonal group O (V) and the main observation is that 
the operator C commutes with the group action. This enables to understand the 
operator C in terms of the representation theory of the orthogonal group and, 
consequently, to compute its spectrum and to describe the associated eigenspaces. 
In this regard, the main technical result of this paper is 

Theorem 3. The operator C admits a discrete spectrum Ai, A2, A„, .., \oo = 
such that 

= (-1)"-^ 
n (n + 1) 

Moreover, dim7i (A n ) = 2n+ 1. 

An immediate consequence of Theorem [3] is that the maximal eigenvalue of C is 
Amax = 1/2, its multiplicity is equal 3 and there exists a spectral gap of Ai — A3 = 
5/12, which separates it from the rest of the spectrum. Now, consider the vector 
space V = TL (A max ) and define the maps 
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The second main result of this paper asserts that the vector space V equipped 
with the collection of tautological embeddings {i x : x G X} is isomorphic to the 
vector space V equipped with the mappings {^iiGl} and, in addition, this 
isomorphism is proportional to the canonical morphism a can :V^>TL which sends 
a vector v G V to the section a can (v) G TL whose value at the point x is the 
orthogonal projection of i> on the plane P^. All of this is summarized in the following 
theorem: 

Theorem 4. The morphism r = y/3/2-a can maps the vector space V isometrically 
onto VcH. Moreover, 

T0i x = ip x , 

for every x G X . 

0.5.1. Proof of Theorems^ and\^ The proof is based on the following three asser- 
tions: 

Assertion 1: The vector space Hn approximates the vector space H. 
Assertion 2: The operator Cjv approximates the integral operator C. 
Assertion 3: The vector space approximates the vector space V. 

The validity of the first two assertions depends on our principal assumption 
that the points in Xn are chosen independently and uniformly at random, which, 
implies that the Haar measure on X is approximated by the (normalized) counting 
measure on Xn- The validity of the third assertion depends also on the existence 
of a spectral gap for the operator C which implies that the maximal eigenspacc 
can be computed in a numerically stable manner. Consequently, Theorems [1] and 
[2] follow from Theorem g] 

Remark 2. The reconstruction of the orthogonal maps {i x : x G Xn} is a non- 
linear problem because of the orthogonality constraint. One of the appealing prop- 
erties of the intrinsic reconstitution algorithm is that it reduces this problem to a 
problem in linear algebra - the computation of the maximal eigenspace of a linear 
operator. Another appealing property is its remarkable stability to noise which, us- 
ing a bit of random matrix theory arguments (see ^\), follows from the spectral 
gap property. Other existing reconstruction methods, like the angular reconstitution 
method (see [9] and ^\), do not enjoy this important stability property. 

0.6. Structure of the paper. The paper consists of three sections except of the 
introduction. 

• In Section [1] we begin by introducing the basic analytic set-up which un- 
derlies cryo-EM. Then, we proceed to formulate the main results of this 
paper, which are: A complete description of the spectral properties of the 
common lines operator C (Theorem [5]) and the admissibility of the intrinsic 
reconstitution algorithm (Theorem [5]) . 

• In Section [2 we prove Theorem [5l in particular, we develop all the repre- 
sentation theoretic machinery which is needed for the proof. 

• Finally, in Appendix \^ we give the proofs of all technical statements which 
appeared in the previous sections. 

Acknowledgement: The first author would like to thank Joseph Bernstein for 
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would like to thank the MPI institute at Bonn where several parts of this work were 
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1. Preliminaries and main results 

1.1. Set up. Let (V, B) be a three dimensional Euclidian vector space over R. The 
reader can take V = R 3 equipped with the standard inner product B st d : R 3 x R 3 -> 
R. Let O (V) = O (V, B) denote the group of orthogonal transformations with re- 
spect to the inner product B; let SO (V) C O (V) denote the subgroup of orthogonal 
transformation which have determinant one; let 9 £ O (V) denote the element —Id. 
Let S (V) denote the unit sphere in V, that is, S (V) = {v £ V : B (v, v) = 1}. 

1.2. The vector bundle of planes. Let — > S(V) be the real vector bundle 
with fibers $j\ x = x 1 - and let TL =T (S (V) ,Sj) denote the space of smooth global 
sections. The vector bundle 9) admits a fiberwise Euclidian structure induced from 
the one on V, which in turns yields a (pre) Euclidian structure on TL (here, the 
word pre just means that TL is not complete). In general, in this paper we will not 
distinguish between an Euclidian/Hcrmitian vector space and its completion and 
the correct choice between the two will be clear from the context. 

In addition, admits a natural O (V) equivariant structure which induces an 
orthogonal action of O {V) on the space of global sections TL which sends a section 
s £ TL to a section g ■ s, given by 

(g ■ s) (x) = gs (g~ 1 x) , 

for every x £ S(V). This makes TL into an Euclidian representation of O (V). We 
will also consider the complexified vector bundle Cf) and its space of global sections 
CTL =r (S (V) ,Cf)). The vector bundle O} is equipped with an Hermitian inner 
product induced from the Hermitian product (— , — ) on CV which is given by 

(u, v) = B (u, v) , 

where (— ) : CV — > CV is the Galois conjugation. Consequently, CTL is a (pre) 
Hemitian vector space supporting a (real) unitary representation of the group O (V). 

1.3. The operator of common lines. We define an integral operator C : TL — > TL 
capturing a basic relation in three dimensional Euclidian geometry. 

The operator C is defined as follows: For every pair of points x,y £ S (V) such 
that x ^ ±y, consider the intersection x 1 - n y 1 - of the corresponding orthogonal 
planes, which is a line in V. There are two tautological embeddings 

C x .y : i 1 !!!/ 1 ^! 1 , 

c y , x : ^n/M/. 

Using these embeddings we can define a rank one operator C (x, y) : Sj\ y — > $}\ x , 
given by the composition C x . y oC y x . The collection {C (x, y) : x ^ ±y} yields a well 
defined smooth section of f) S F>* on the complement of the union of the diagonal 
and the anti-diagonal. It is not difficult to verify that this section extends to a 
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distribution section C e V (S (V) x S (V) ,$) Klf)*) which, in turns, establishes a 
kernel for an integral operator C : Ti — > Ti given by 

C(s)(x)= J C(x,y)(s(y))dy, 
yeS(V) 

for every s G Ti, where we take dy to be the normalized Haar measure on the sphere. 
Since C (x, y) — C (y, x) f , this implies that C is a symmetric operator. In addition, 
it is evident that C commutes with the O (V) action, namely C (g ■ s) = g ■ C (s) 
for every s &Ti and g S O (V) . 

The operator C is referred to as the operator of common lines and the main 
technical part of this paper will be devoted to the investigation of this operator. 

1.4. Main results. Our goal is to describe an intrinsic model of the Euclidian 
vector space V which is expressed in the terms of the common lines operator and 
the Euclidian structure of Ti alone. 

The main technical result of this paper is 

Theorem 5. The operator C admits a discrete spectrum Ai, A2, .., A„, ..Aoo = 0, 
such that 

(-If 1 
An ~ n(n+lY 
Moreover, dim7i (A„) = 2n+ 1. 

For a proof, see Section [5) 

Intrinsic model: Take V = Ti (A max ) to be the maximal eigenspace of C. 
There are two immediate implications of Theorem [5] that we will consider: 

• The vector space V is three dimensional. 

• There exists a spectral gap of Ai — A3 = 5/12 which separates A max from 
the rest of the spectrum. 

Fix r = Let ev x : Ti — > x 1 - denote the evaluation morphism at the point 

ieS (V). For every x € 5* (V), define the morphism 

ip x = r~ x ■ (ev x \y) :x ± ^W. 

Let a can '■ V — > Ti be the canonical morphism which sends a vector v G V to the 
section a can {v) S Ti, defined by 

Ctcan (v) (x) = Pr x (v) , 

where Pr^ is the orthogonal projection on the plane x . The morphism a can is a 
morphism of representations of the orthogonal group O (V). 
Finally, define 

t = r • a can :V^Ti. 
The morphism r identifies the Euclidian vector space V together with the tauto- 
logical embeddings {i x : x 1 - — * V} with the Euclidian vector space V together with 
the maps {ip x : x 1 - — > V}. All of this is summarized in the following theorem: 

Theorem 6. The map t maps V isometrically onto V. Moreover, 

T0i x = Lp x , 

for every x 6 S (V). 
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For a proof, see Appendix [A] (the proof uses the results and terminology of 
Section El). 

2. Spectral analysis of the operator of common lines 
2.1. Set-up. It will be convenient to extend the set-up a bit. 

2.1.1. Auxiliary vector bundles. We introduce the following auxiliary vector bundles 
on S (V). Let 91 — > S(V) be the vector bundle of normal lines with fibers yi\ x = Mx 
and let JV = T (S (V) , 91) denote the corresponding space of global sections. Let 
Vs(v) denote the trivial vector bundle with fiber at each point equal V and let V 
= T (S (V) , V s{v) ) = T®V where T = C°° {S {V) , R), we have 

V = TL®N. 

All the vector bundles are equipped with a fiberwise Euclidian structure which 
is induced from the one on V and consequently the spaces of global sections are 
Euclidian. In addition, all the vector bundles are equipped with a natural O (V) 
equivariant structure which is compatible with the Euclidian structure and conse- 
quently the spaces of global sections form Euclidian representations of the group 
0(V). We will consider these spaces as representations of the subgroup SO (V) 
and remember also the action of the special element 9 £ O (V) which commutes 
with the action of SO (V). 

2.1.2. The operator of orthographic lines. We define an integral operator O : TL — > 
TL which we refer to as the operator of orthographic lines. This operator captures 
another basic relation in three dimensional Euclidian geometry which, in some 
sense, stands in duality with the common lines relation. 

The operator O is defined by the following kernel: For every pair of points 
(V) such that x ^ ±y, consider the pair of unit vectors 

„ _ Pr.(y) e , 



|Pr* (y)ll 



Py v (x) ^ 



HPr*(aO| 

In words, the vector o XjV is the normalized projection of the vector y on the plane 
x and similarly the vector o y x is the normalized projection of the vector x on the 
plane y . Define a rank one operator O (x, y) : y — > x by 

O (x, y) (v) = B (o ViX ,v) o XtV , 

for every » 6 j 1 . The collection {O (x,y) : x ^ ±y} yields a well defined smooth 
section of S) H Sj* on the complement of the union of the diagonal and the anti- 
diagonal which extends to a distribution section O € T' (S (V) x S (V) , Sj M S)*) 
which, in turns, yields a symmetric integral operator O : TL — » TL which commutes 
with the O (V) action. 

Given a pair of unit vectors x ^ ±y, the following observations are in order. 

• The orthographic lines M.o x , y £ x and M.o y . x £ y are orthogonal to the 
common line x D y . 

• The kernel 0(x,y) satisfy 0(x,—y) = 0(—x,y) = —1 • 0(x,y) which 
means that O {x,y) depends on the choice of the unit vectors x,y and not 
only on the planes x ± ,y ± . This should be contrasted with the analogue 
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property of the kernel C(x,y) which satisfies C(x,—y) = C (—x,y) = 
C(x,y). 

2.1.3. The operator of parallel translation. We define the integral operator T = 
C — O : H — > TL. For every x,y £ S (V), such that x ^ ±y, the kernel T (x, y) : 
y — > x is a full rank operator and it is not difficult to verify that T (x, y) is the 
operator of parallel translation along the unique geodesic (large circle) connecting 
the point y with x. Consequently, we will refer to T as the operator of parallel 
translations. 

The strategy that we are going to follow is to study the spectral properties of 
the operator T , from which, as it turns out, the spectral properties of the operators 
C and O can be derived. 

2.2. Isotypic decompositions. The spaces TL, V,Af and T form Euclidian repre- 
sentations of the group SO (V) and as such decompose into isotypic components^ 

oo 

n = 0w„, 

n=0 
oo 

n=0 
oo 

V = 0K, 

n=0 
oo 

J~ = J~ n > 

n=0 

where we use the subscript n to denote the isotypic component which consists of 
the unique irreducible representation of SO (V) of dimension 2n + 1. In addition, 
the element 6 £ O (V) acts on all these spaces, thus decompose them into a direct 
sum of two components 

n = n + ®n-, 

M = J\f + © N~ , 
V = v+©v~, 

where we use the superscript + to denote the component on which 9 acts as Id and 
the superscript — to denote the component on which 9 acts as —Id. We will refer 
to the + component as the symmetric component and to the — component as the 
anti-symmetric component. 

The following theorem summarizes the properties of these decompositions which 
we will require in the sequel. 

Theorem 7. The following properties hold 

(1) Each isotypic component T n is an irreducible representation. Moreover, 
T n = when n is even and T n = T~ when n is odd. 



We remind the reader that an isotypic component is a representation which is a direct sum 
of copies a single irreducible representation. 
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(2) Each isotypic component Af n is an irreducible representation. Moreover, 
■Mi = Afn when n is even and Af n = N~ when n is odd. 

(3) The isotypic component Ti$ = and each isotypic component 7i n , n > 1 
decomposes under 9 into a direct sum of two irreducible representations 
Tin © Ti-n ■ 

(4) The isotypic component Vq is equal to the symmetric trivial representation 
1 + and each isotypic component V n , n > 1 decomposes under 9 into a direct 
sum of three irreducible representations Ti.^ © Ti.~ © Af n where 

„ _ J + n even 
\ — n odd 

For a proof, see Appendix 

Since the operators C, O and T commute with the O (V) action, they preserve 
all the above decompositions. 

Proposition 1. The following properties hold: 

• The operator T acts as scalar operator on TL n , moreover, T\f^ n = X n Id 
where X n ^ 0. 

• The isotypic component TC^ C ker C , moreover, Cy^- — X n Id . 

• The isotypic component Tl~ C kerO, moreover, 0| W + = —\ n Id. 

For a proof, see Appendix [A] 

The rest of this section is devoted to the computation of the eigenvalues A n . 

2.3. Computation of the eigenvalues. Fix a point x G S (V). Let T x — 
{g G SO (V) : gx = x} be the subgroup of rotations around x. Let (H,E,F) G 
CLie (SO (V)) be an an s^ triple associated with T x . 

2.3.1. Spherical decomposition. For each n > 1, the complexified (Hilbert) space 
C7i rl admits an isotypic decomposition with respect to the action of T x 

n 

cn n = 7C, 

m=—n 

where H acts on "H™ by 2mld. Since CH n = C7i+ © CH~, each H™ is two 
dimensional. 

2.3.2. Strategy. Given a section u n G CH n , by Proposition [1] Tu n = \ n u n . If, in 
addition, u n G W„ then, as will be shown, u n can be chosen such that u n (x) ^ 0. 
Under this choice, the eigenvalue A„ can be computed from the equation 

A„ (u n (x) , u n (x)) = (u n (x) , Tu n (x)) . 

The following proposition gives an explicit formula for (u n (x) ,Tu n (x)). But 
first, we need to introduce one additional terminology. 

Fix a unit vector l G S (x J C V and let A[ G Lie (T; ) be a vector such that 
the morphism exp : [0, 2n) — > Ti given by exp (9) = e eAl o is an isomorphism. 

Proposition 2. The following equation holds: 

7T 

(2.1) A„ (u n (x) , u„ (x)) = /" /x (6) (u n (x) , e- Mi o w „ ( e 9A 'o x ))d9, 



where ji (9) = sin (9) jl. 
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For a proof, see Appendix [A] 

Our strategy is to construct "good" section u n £ H n and then to use Equation 

(E3). 

2.3.3. Construction of a "good" section. For every n > 0, choose a highest weight 
vector ip n £ CJF„ with respect to (H, E, F) ( Hip n — 2mp n ). In addition, choose a 
highest weight vector vi £ CV {Hv\ = 2vi). 

For every n > 1, the section ip n _i ® f 1 is a highest weight vector in CV„. In 
order to get a section in TL nl first step, apply the lowering operator F and consider 
the section u n = F^ 1 (V> n _i O «i) £ V, 1 ,. 

Let us denote by P„_i the weight zero spherical function F n ~ 1 ip n _ 1 £ and 
note that, under appropriate choice of coordinates, P n ~i is the classical spherical 
function on the sphere S (V). The following proposition gives an explicit expression 
of u n in terms of the function P n -\ and the vector v\. 

Proposition 3. The section u n can be written as 

(2.2) u n = P„_! ® vi + -EP n -i (8) Fvi + —^——E 2 P n ^ 1 ® F 2 Vl . 

n 2n(n + 1 ) 

For a proof, see Appendix [AJ 

The second step in order to get a section in 7i n is to define 

Un (y) = PyU n (y) , 

for every y E S (V), where P y is the orthogonal projector on y- 1 . 

Proposition 4. The section u n is symmetric or anti- symmetric depending on the 
parity of n as follows 

u n G T^n' 1 w hen n is even, 
u n £ Hn' 1 when n is odd. 

For a proof, see Appendix |A] 

Now, we are ready to finish the computation. Using Equation (|2.ip , we can write 

7T 

A„ {u n (x) ,u n (x))= J (i (6) (e eA 'o Un ( x ) , u n (e eA 'o x ))d6. 



Considering formula (|2.2[) . it is evident that the functions EP n -\ and E 2 P n -\ 
must vanish at a; £ S (V) since these are function of weight different then zero with 
respect to the action of T x . Hence 

U n (x) = U n (x) = P n -1 (x) Vi £ Cx ,l . 

For k = 0, 1, 2, define the integrals 

ll = {Un[x) \ n{x)) J ,(0)(e^Un( X ),F^(e^ X )F^)dO 



7T 

= o An „2 / ^) Ek Pn-i (e eA '°x) (e 0A 'ov u F k Vl )de, 

P n -1 (X) \\Vi\\ J 
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and note that the eigen value A„ can be expressed as 
(2-3) A„ = J° + hi < ' 



n n 2n(» + l) 

Theorem 8 (Main technical statement). For k = 0, 1, 2, the integrals J„ are egwai 
to 





f 1 n = 1 






1 1/6 n = 2 , 






[ n > 3 






f n=l 




/,! = | 


-2/3 n = 2 






[ n > 3 






f n 


= 1 


^ - 1 


1 n 


= 2 




[ 2(-l)"~ 1 n 


> 3 



For a proof, see Subsection 12.41 

Consequently, using Theorem[8]and Equation ()2.3|1 we obtain the desired formula 

(_!)»-! 

An — 



1) 

which proves Theorem [5J 



2.4. Proof of the main technical statement. Let (ei, e^, ez) be an orthonormal 
basis of V. Silently, the reader should think of the basis vector e$ as standing in 
place of the fixed unit vector igS (V) and of the basis vector ei as standing in 
place of the vector Iq e S'(x ± ). It is possible to choose vectors A ei 6 Lie(T ei ) 
which satisfy the relations 



[A e3 , A ei ] 


— A 

^e 2 i 


[A e3 ,A e . 2 ] 


= ~Ae 


[A ei , A e . 2 ] 


— A e3 , 



and, in addition, satisfy [A ej ,Ae 3 .] = A Bi ej for every 1 < i,j < 3. We can define 
the following triple (H, E, F) which is associated with T e3 



H = -2iA e3 , 
E = iA e2 -A ei , 
F = A e± -{-iA e2 . 

2.4.1. Spherical coordinates. We introduce spherical coordinates / : (0, 2tt) x (0, w) 
S (V) given by / (ip, 9) = g v ■ (cos (9) e% + sin (9) ei) where 

(cos (tp) — sin (</?) 
sin (tp) cos (tp) 
1 



REPRESENTATION THEORETIC PATTERNS IN THREE DIMENSIONAL CRYO-EM 13 



In the coordinates (<p, 9) the operators H, E, F are given by the following formu- 
las 

H = 2id v , 

E = -e~ llp {idg + cot (0) d v ) , 
F = -e iv (idg - cot (0) d v ) . 

2.4.2. Highest weight vector in V . The vector v\ = —e% + ie2 is a highest weight 
vector in V and we note that \\vi || 2 = 2. For k = 0, 1, 2, let us denote by j k (9) the 
function (e e »i>i, F k v\j. Explicit calculation reveals that 

j°(9) = cob(0) + 1, 
j l {8) = 2ism(9), 
j 2 (9) = 2 cos (6)- 2. 

2.4.3. Spherical function in T n . For n > 0, let P n £ ^ denote the unique weight 
zero spherical function which satisfies the normalization condition P n {e^) = 1. 
Define the generating function 

oo 

n=a 

The generating function G admits an explicit formula. 
Theorem 9 ([8]). 

(2.4) G(ip,9,t) = (l -2t cos (9) +t 2 )~ 1/2 . 

-l °° 

Remark 3. Note that G (0, 0,i) = (1 — t) — ^2 t n which is compatible with the 

n=0 

normalization condition P n (0,0) = 1. 

Applying the raising operator E we obtain the generating functions 

oo 

EG(<p,B,t) = £EP„(M)* n , 

oo 

E 2 G(^,9,t) = ^£ 2 P„(^,0)i". 

n=0 

Granting formula (|2.4p . explicit calculation reveals that 

EG(ip,6,t) = ie- lv t sin {9) (l-2t cos (9) + t 2 y 3/2 , 
E 2 G(<p,9,t) = -3e- 2lv t 2 sin 2 (9) (l - 2t cos (9) +t 2 )~ 5/2 . 

2.4.4. Putting everything together. In terms of our choice of the highest weight 
vector v\ and the spherical functions -P„'s, the integrals /„, k = 0, 1, 2, are given in 
the spherical coordinates (</?, 9) by 

7T 

J* - 1 / p(6)E k P n - 1 (p,0)j k (0)dO 

Pn-l{x)\\vi\\ J 

7T 

= i J f x(6)E k P n _ 1 (0,9)j k (6)d6. 
o 
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For k = 0, 1, 2, define the generating functions 

n=Q 

Each I k (t) can be expressed as the integral 

7T 

I h (*) = ^ / E k G(0,6,t)j k (6)d6. 
o 

Explicit calculation of the integrals I fc (t) reveals that 

!•(«) - i(i + 5< 

"< f > - i (4) ■ 

1 OO 

J 2 (i) = -(4(l + t)- 1 -4t-4)=2^(-irt 

n=2 

From the above formulas, using Equation (|2.3|) .we get 

( _ ir -i 



n (n + 1) 

for every n > 1. This finish the proof of Theorem [5J 

Appendix A. Proofs 

A.l. Proof of Theorem [6J Since r is a morphism of O (V) representations, it 
maps V isomctrically onto TC^ - the unique anti-symmetric copy ( 9 acts by —1) of 
the three dimensional representation of SO (V), which, by Proposition [TJ coincides 
with V = H(A max ). 

Evidently, r is an isometry, up to a scalar. Hence, it is enough to show that 
Tr (t o r T ) = 3, which we verify as follows: 

Tr (r o r T ) = r 2 - Tr (a can o a^ an ) = - J Tr (i£ o i x ) dx 

xes{v) 

3 
2 

x£S(V) 

Finally, the relation t o i x = ip x follows from 

eX x \Y ° Ctccm = P? x - 

This concludes the proof of the theorem. 

A. 2. Proof of Theorem[7l Property 1 is the classical result of spherical harmonics 
on the two dimensional sphere, which can be found for example in [SJ . We just 
note that the representation T n consists of the restriction to S (V) of harmonic 
polynomials of degree n, which implies that T n = when n is even and T n = T~ 
when n is odd. 

Property 2 follows from property 1 since Af can be trivialized using the O (V) 
invariant section s S Af where s (y) = y, for every y 6 S (V). 



2dx = 3. 
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We now prove Properties 3 and 4 simultaneously. 

Since V = T © V as a representation of 0(V), we can compute the isotypic 
components of V in terms of the isotypic components of T . The computation 
proceeds as follows: 

• For n = 0, JF ® V = V~. 

• For n > 1, T n ® V = © 7)*_ x © (JF„ © 7)* © (JF„ © 7)' l+1 where 

? j + n odd 
1 — n even 

The decomposition of T n © 7 as a representation of SO (7) is computed 
using the branching rules of a tensor product and the action of 9 € O (7) 
is derived from the facts that V = V~ and Property 1. 

This implies that the isotypic components of V are 

• For n = 0, V„ = 1+. 

• For odd n > 1, V» = {F n -i © V)~ © (J>, © 7)+ © (.Fn+i © V)~. 

• For even n > 1, V„ = © 7)+ © (JF„ © V)~ © © 7)+ 

Combined with Property 2 and the fact that V = Ti © A/" yields Properties 3,4. 
This concludes the proof of the Theorem. 

A.3. Proof of Proposition [TJ Fix n > 1. Denote ft = H n and ft± = ft±. The 

statement that TL + C kerC follows from the facts that C (ac, — y) = C (—x,y) — 
C (x,y) and that a section s e H + satisfies that s (— x) — —9(s)(x) — —s(x). 
Similarly, the statement that H~ C kerO follows from the facts that O (x, —y) — 
O (—x, y) — —O (x, y) and that a section s S ft - satisfies that s (—x) = —9 (s) (x) = 
s{x). 

Since, by definition, T = C — O, this implies that Cyn- — T\h- and — Oi^+ = 
T\ji+. Since T commutes with the action of SO (7) and TL are irreducible repre- 
sentations 

T\ n ± = X ± Id. 

We are left to show that A + = A© The argument proceeds as follows: 

Let us denote by ev x : Cft — > Cx 1 ' the evaluation map at the point x. Since x 

is fixed by the group T Xl this implies that ev x is a morphism of representations of 

T x . Moreover, ev x induces an isomorphism of weight spaces 

ev x : H ±A -=-> ar 1 ' 1 . 

Fix sections u G TC^' 1 normalized to be of norm 1. Applying Formula (|2.ip . we 
can write 

7T 

(A.l) A ± (ew 2; M ± ,e^u ± ) = J n (9) (ev^, ev x (e M, »u ± )) d6 

o 

= (ew a; u ± , ev x (7r ± (p) it n )) , 

where tt ± : T x -> U (CH*) are the group actions restricted to the subgroup T x and 
/I is the function on T x corresponding to /i via the isomorphism e 'o . 
Equation (|A. 1|) . implies that 

A ± = (u ± ,ir :t (m)m ± ) cw± . 
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This implies that A are characterized solely in terms of the irreducible repre- 
sentation tt^~ : SO (V) — » U (CH ± ), which, in turns, implies that A + = A~. 
This concludes the proof of the proposition. 

A. 4. Proof of Proposition [2j Let / : T x x (0, n) —> S (V) be the spherical coor- 
dinates on S (V) given by / (g, 9) — ge eA 'ox. In these coordinates, the normalized 
Haar measure on S (V) is given by dg El fx (9) d9 where dg is the normalized Haar 
measure on T x and \x (9) = sin (6) /2. 

The section u„e TL\ is a character vector with respect to the group T x , let us 
denote this character by x '■ T x — ► S and note that we have g ■ u n = x (d) u n , for 
every g G T x . Now, compute 



A„ (un (x) , u n (x)) = (u n (x) , Tu n (x)) 

(u n (x),T(x,y) u n (y))dy 



yes(V) 

„8A ln \„. I OA 



dgjfi {9) (u n (x) , T (x, ge 9Al o x ) Un (ge eA >o x) ) d9 
o 

7T 

dg J n(d) (u n {x),gT (x, e eA 'o x) g- l u n {ge eAl «x))d9 
o 

dg J n (9) (g-'un (x) , T (x, e 9Al « x) g~ x u n (ge eAl ° x) ) d9 
o 

7T 

dg [ p (9) (u n (x) ,T(x, e eAl <> x) u n (e eAl ° x) ) d6 



li{9) (u n (x),e- 0Al o Un (e eA '°x))de, 

where, step 4 follows from the fact that T commutes with the action of SO (V) which 
is equivalent to the property that T (gx, gy) = gT (x, y) f? -1 , for every x,y G S (V) 
and g € SO (V) which, in particular, implies that T (x, ge eAl o x) =T {gx, ge eAl ° x) = 
gT (x, e eAl ox) g~~ x and step 7 follows from the fact that T (x, e eAl ox) is the opera- 
tor of parallel translation along the big circle connecting the point e eAl o x with the 
point x. 

This concludes the proof of the proposition. 

A. 5. Proof of Proposition [3j First we note the following simple fact: The op- 
erator EF : CV n — > CV n preserve the weight spaces V„, and, moreover, since CV„ 
is a representation of highest weight 2n with respect to the sl% triple (H, E, F) we 
have 

(A.2) EF \Vi = (n + l)(n-l + i) Id, 

for I = —n, ... n. 
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Now, calculate 

Un = F- 1 ® vi) = £ (" 7 ^ ^""^ ® (V»„-l ® «l) ■ 

i=0 ^ ' 

Since is a representation of highest weight 2 with respect to the s?2 triple 
(H, E, F), all tensors of the form (— ) C3> P i>i, for fc > 3, vanish. This implies that 
the above sum is equal to 

F^V„-i ® «i + (n - 1) P-Vi ® Fv x + ( n - 1 )( n - 2 ) f n-3^_ i g 

Recall that P n -i = F n ~ x ij) n _ l . Explicit calculation, using formula (|A.2[) . reveals 
that 

P n ~V„-i = p 1 rr SPn-l, 

P™ 3- 0n-i = 7 — 1 — TT^E 2 ^- 1 - 

(n — 2) (n — 1) n {n + 1) 

Combining all the above yields the desired formula for u n . 
This concludes the proof of the proposition. 

A. 6. Proof of Proposition^ The statement follows directly from the facts that 
V = V~ which implies that 9 (F k vi) = —F k v\ and that P n -i S T' n _ x where 



? 



f n odd 
- n even 

This concludes the proof of the proposition. 
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